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Abstract 

The asymptotic expansion of the massive scalar field propagator on a n- dimensional 
lattice is derived. The method used is based on the evaluation of the asymptotic 
expansion of the modified Bessel function I v (y 2 f3) as the order v grows to infinity. 
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1 Introduction 



In perturbation theory, the advantage of a coordinate space description as a 
way of studying the divergences arising in Feynman diagrams has been dis- 
cussed by many authors (e.g. [1,2]). In the continuum, the analytic expression 
for the scalar field propagator in position space is well-known. On the lattice, 
however, the standard representation for the scalar propagator involves inte- 
grals over Bessel functions and has proved to be very difficult to analyse in 
the continuum limit. 

Several attempts have been made to derive a suitable expansion for the lattice 
scalar propagator in the limit where the lattice spacing goes to zero. In partic- 
ular, we mention the analysis performed by Luscher and Weisz for the massless 
propagator [3] as well as the research carried out by Burgio, Caracciolo and 
Pelissetto in [4] . In the first of these studies, the x- dependence of the massless 
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propagator is asymptotically derived. In the second, the m-dependence is ob- 
tained for the propagator at x — 0. The propagator at any other point i^O 
is then expressed, through a set of recursion relations, in terms of the value 
which its massless counterpart assumes on an element of the unit hypercube 
(i. e, Xn either or 1). In the present paper we wish to tackle the most general 
case of both m and x non-vanishing. 

The procedure that we adopt attacks the question directly at its core. We 
derive an asymptotic expansion of the modified Bessel function of the first 
kind which appears in the expression of the lattice propagator. The order and 
argument of this Bessel function go to infinity at different rates as the lattice 
spacing decreases to zero. We have found no tabulated asymptotic expansion 
for this particular case. Expansions for the modified Bessel function of the 
first kind, I u (x), are, indeed, available for the cases where either the argument 
x or the order v becomes large, or for the case where both x and v grow to 
infinity keeping the value of their respective ratio constant (approximantion 
by tangents). Unfortunately, none of the cases just mentioned characterizes 
the modified Bessel function at hand, and we have had to develop, as a result, 
a new asymptotic expansion. 

As a perturbative application of the technique developed in the analysis of 
the continuum limit expansion of the lattice scalar propagator, in the closing 
section of this paper we consider the mass renormalization of the discrete A0 4 
theory in coordinate space. 

We introduce now briefly the notations that will be assumed in the following 
sections. Throughout this paper, we shall work on a n-dimensional lattice of 
finite spacing a. The convention 

n 

x p = x^, with p an integer , (1) 



will also be used. 

The Euclidean free scalar propagators in a n-dimensional configuration space 
will be denoted by A c (x; n) and A L (x; n) with 
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referring to the propagator evaluated in the continuum and on the lattice, 
respectively. Note that in eq. (3) we have introduced the short-hand notation 




2 The propagator on the lattice 



It is a well- known result that the continuum scalar field propagator in a 4- 
dimensional configuration space can be written in terms of the modified Bessel 
function K x . The technique consists in using the Schwinger's representation for 
the propagator in momentum space. Re-expressing, then, the integral over the 
four-momentum as a product of four one-dimensional integrals and completing 
the square in the resulting exponential, we can finally perform the integration 
and obtain our final formula. Adopting the notational convention given by eq. 
(1), we write the n-dimensional expression as follows 
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The derivation of the standard representation for the lattice propagator in 
configuration space is carried out in a much similar fashion. Indeed, we have 



A L (x;n)= f-^- 
K 1 ' J 2tt " 
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with IfM (jpj corresponding to the modified Bessel function of the first kind. 
Unfortunately, the integral appearing in eq. (6) cannot be trivially solved. 
As a consequence, we are not able in this case to express the propagator in 
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closed form. What we really wish to do here, though, is to show that in the 
continuum limit, i.e. when a — > 0, A L (x;n) is given by the sum of A c (x;n) 
plus a series of correction terms depending on increasing powers of the lattice 
spacing a. 

As already mentioned in the introduction, the most direct way to proceed in 
order to achieve our goal is trying to derive the asymptotic expansion for the 
modified Bessel function I*m (%) as a — > 0. 

The strategy to adopt for the actual derivation of the expansion is determined 
by the fact that the order and the argument of the modified Bessel function 
become large at different rates. The standard techniques of global analysis 
(e.g, the steepest descents method) are not of much use in this case. As a 
consequence, we are forced to choose here a local analysis approach. This 
implies beginning our study by examining the differential equation satisfied 
by the modified Bessel function at hand. With the purpose of determining 
uniquely the asymptotic behaviour of the solution, we shall also impose in the 
end the condition that the series representation derived reproduce (through 
its leading term) the continuum result A c (:r;n). 
We first commence by setting 

P = and v = — , (7) 

xi a 

H 1 



i.e. l xJa (2a/a 2 )^l u (u 2 p). 

Thus, we wish to find an expansion for \ u (v 2 P) as v — > oo. We observe now that 
\ v {y 2 P) satisfies a differential equation of the form: x 2 ^+x^ — (x 2 + u 2 )y = 0. 
Hence, performing the change of variable 



x 



v 2 P 
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we obtain 
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I„(v 2 P) = 0. 



(9) 



Eq. (9) simplifies slightly if we make the substitution 



(10) 



C represents here the free parameter whose value shall be fixed later on ac- 
cording to the prescription made at the beginning of this paragraph. 
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By using eq. (10), we now get a new differential equation for Y(f3, v), namely 
d 2 Y(P, v) 



dp 2 



3 2 



Y0M = o. 



(11) 



We consider at this stage the limit v — > oo. We now need to make some 
kind of assumption on the form of leading term governing the expansion of 
the solution to eq. (11) as v becomes large. With this purpose, we consider 
a substitution originally suggested by Carlini (1817), Liouville (1837) and 
Green (1837) and whose underlying idea is that the controlling factor in the 
asymptotic expansion of a function is usually in the form of an exponential. 
Hence, we assume the asymptotic relation 

Y((3,u) ^ e » 2 P+jMMd^ ( 12 ) 



Actually, the form which we have used here for Y(j3, v) differs slightly from the 
standard Carlini's substitution since it implicitly assumes that not only the 
leading term, but also all the subsequent terms in the expansion are expressible 
as exponentials. The reason which led us to adopt eq. (12) as a viable option 
is essentially twofold: on one hand, it is dictated by the necessity of recovering 
the well-known expression for the free scalar propagator in the continuum; 
on the other, the idea spurred in analogy to the method adopted in Watson 
[5] and originally due to Meissel (1892) in order to deduce the asymptotic 
expansion for the Bessel function of the first kind, 3 v {yz) with v large. 
Carrying out the substitution in the differential equation for Y(j3, u), we end 
up with a first-order inhomogeneous differential equation for the new function 
w(j3, u), i.e. 

9 -^l + w 2 ((3, v) + 2^(/3, v) - ^-1 ~ 0. (13) 



Our aim, at this point, is to look for a suitable series representation for the 
solution to this equation. With this purpose, we make the following observa- 
tions. First of all, if we suppose that 2v 2 w(f3, v) — v 2 / f3 2 ~ represents a good 
approximation to eq. (13) when v — > oo, then it is straightforward to derive 
the leading behaviour of w(j3, v) as {2f3 2 )~ l . We note that this gives exactly 
the continuum result A c (x;n) once the overall solution to eq. (9) is normal- 
ized by fixing the arbitrary parameter C to be equal to {2nu 2 )^ 1 ^ 2 . Secondly, 
we observe that all the non-leading terms in the expansion of w((3, v) should 
feature only even powers of z/, since only even powers of this variable appear 
in eq. (13). 

As a result of the remarks just made, we assume that w(j3, v) admits an 
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asymptotic series representation for v large of the type 

oo 



-2n 



(14) 



n=0 



This is consistent with the approximation that we considered when we deduced 
the form of the leading term of the solution. Substituting now this formula 
into eq. (13), we get 



oo 

2 £ a n+1 (py 

n=-l 



V 2 - i 

-2n ^ V 4 



P 2 



n=0 



J2 a m (P)a n -m(P) 

.m=0 



Matching, order by order, the coefficients corresponding to the same power 
of v, we are now able to deduce the set of relations defining the coefficients 
a n ((3). We have 
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The coefficients a n (/3) can be, at this stage, computed iteratively leading to 
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The asymptotic expansion of l u (v 2 {3) as v — > oo can be now easily derived. 
The result reads as follows 
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Dividing this equation by \ u (v 2 f3) and plotting the result against (5 for a fixed 
(large) value of u, we see a fluctuation of the ratio around unity. This repro- 
duces exactly the behaviour which ought to be expected by the ratio of two 
functions asymptotic to each other. Furthermore, we observe that the fluctu- 
ation around 1 is a feature of the whole positive /3-axis, therefore implicitly 
suggesting the validity of our expansion even for large values of f3. 
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Fig. 1. Behaviour of the ratio of the asymptotic expansion of l v (v 2 (3) with l v {u 2 (5) 
as a function of with v equal to 100. 

Recalling the expressions of (3 and v in terms of a, x M and a, we are now able 
to derive the formula for the n- dimensional propagator up to sixth order in 
the lattice spacing. In fact, we can do more than that. We consider at this 
point a generalisation of A L (x;n) by introducing an arbitrary exponent q in 
the denominator of eq. (6); that is, we now look at 
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The main reason for considering this generalised quantity rests in the fact that 
eq. (20) represents a key element in the expression of general one-loop lattice 
integrals with bosonic propagators and zero external momenta [4]. 
Having reached a formula for the asymptotic expansion of the modified Bessel 
function as the lattice spacing vanishes, the study of the continuum limit of 
A L (x; n; q) is not technically more difficult than the analysis of the same limit 
for the n- dimensional lattice propagator A L (x;n) = A L (x;n; 1). Indeed, we 
simply have now to substitute eq. (19) into eq. (20) and carry out the product 
over the dimensional index /i. The resulting a-integrals are all well-defined 
and finite. We can, therefore, proceed to their evaluation and obtain 



A L (x;n; q) ~„^ 
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The full expression of each of the coefficients A^(x; n; q) (i = 0, 1, 2, 3) is given 
in terms of the new function P p (m; x) defined as 
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representing, as usual, the modified Bessel function of the 



second kind and p a real number. Therefore, we have 
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The expansion obtained for A L (x;n;q) clearly shows how the finite correc- 
tions introduced by formulating the theory on a lattice can be analytically 
expressed by a series of increasing (even) powers of the lattice spacing a with 
coefficients given by analytic functions of the mass and space coordinate times 
a modified Bessel function of the second kind of increasing order p. 
We intend now to demonstrate how eq. (21) is in perfect agreement with both 
the studies performed in [3] and [4]. With this aim, we analyse A L (x;n;g) 
in the limit m(x 2 ) 1 / 2 — > 0. Given the functional dependence of eq. (21) on 
P p (m;x) and given the definition in eq. (22), this translates into considering 
the appropriate expansion for the Bessel function K p m(x 2 )^ . We wish to re- 
call at this point that the series representation of the modified Bessel function 
of the second kind assumes different forms depending on whether the order p 
is a real (p re ) or integer (pm)- In particular, for mix 2 ) 1 ! 2 small and p = p TC we 
have 
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while, for m(x 2 ) 1//2 still small and p = p m (p m ^ 0) Q we find 
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with ip denoting the ip- function [6]. 

Using the relation T(l — p re )r(p re ) = it/ sinp re 7r, we find that P Pro (m;a;) and 
P Pin (m; x) assume the same functional form 2 p ~ 1 T(p)/(x 2 ) p (p > 0) as m — > 0. 
The limits m — * and p rc — > p in are, therefore, uniform. After performing the 
final mechanical analysis to set n = 4 and q = 1, we find that the massless limit 
of our continuum expansion for A L (x; n; q) reproduces exactly the expansion 
obtained by Luscher and Weisz for the massless 4-dimensional propagator [3] 
The limit x — ► proves to be more difficult to analyse. The short- distance 
behaviour is, indeed, singular and the limit not uniform in this case. Note 
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that this observation matches the analogous remark made in [1] about the 
behaviour of the dimensionally regularised propagator. Observe also that the 
logarithmic mass-behaviour described in [4] by Burgio, Caracciolo and Pelis- 
setto is recovered, in our formulation, for integral values of p. 



3 The tadpole diagram in A0 4 theory 



As a direct implementation of the result obtained for the continuum expansion 
of the scalar propagator, we now wish to derive the one-loop renormalization 
mass counterterm of A0 4 theory through the study of the lattice tadpole dia- 
gram. 




Fig. 2. Self-energy tadpole diagram in A0 4 theory. 

In a n-dimensional continuum space, the contribution to the full propagator 
associated with the graph in Fig. 2 is well-known and proportional to the inte- 
gral / d n zA c (x — z)A c (0)A c (z — y). The ultraviolet behaviour of the diagram 
is entirely due to the divergence in A c (0) [1]. 

In terms of our notational conventions, the lattice version of the tadpole graph 
in n-dimensions is immediately written as 
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tad 



71 ) 



(29) 



with A the coupling constant and fi^ n a multiplicative dimensional factor 
introduced to preserve the dimensional correctness of the theory. Note that, 
in eq. (29), the lattice propagators have replaced their continuous counterparts 
and the summation over the lattice sites has taken the place of the continuum 
integral. 

Our present goal is to examine eq. (29) as a — > 0. Using the asymptotic 
expansion of the propagator as derived in eq. (21), it is straightforward to 
obtain up to fourth order in the lattice spacing 
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a 4 



h(x,y,z)A^O;n)+g(x,y,z)A^(0;n) + f(x,y,z)A L 4 (0;n)]+...} , (30) 



with f(x, y, z), g(x, y, z) and h(x, y, z) denoting functions which are associated 
with products of the type A\ Aq, A2 Aq and A\A\ plus AJ; A^, respectively. 
We observe that the n- dimensional coefficients which appear in the expressions 
of /, g and h are exclusively evaluated at x—z and z—y and correspond, hence, 
to infinities which are, ultimately, integrable. As a result, the singularities 
in M|? ad are fundamentally generated by the poles in Aq(0;ti), A^Ojn) and 
A 4 (0;n) only It is now of paramount importance to remark that each of 
the latter quantities scales, implicitly, with the lattice spacing a. This scaling 
needs to be made explicit in the analysis by considering the mass as physical, 
i.e. mR = ma, and by recalling that, in a discrete formulation, the space 
coordinate is also defined in terms of the lattice spacing. 
At this stage, the isolation of the UV-divergences can take place along the 
lines of the method introduced in Collins for the study of the continuum 
case. Indeed, performing a point- splitting of the tadpole z-vertex through the 
introduction of an arbitrary variable £r = e/a, such that e R — > (a fixed), 
the investigation of the divergent behaviour in each of the A^Ojn) (i = 
0, 1, 2) translates now into examining the divergences in A^Er; m^; n) for 
vanishing e R and n fixed and non-integral. We note that each of the series 
expansions obtained splits, naturally, into two sub-series, with the poles in 
£r all contained in the first. In terms of e R , the second sub-series is, in fact, 
analytic. Thus, focusing on the singular contribution, we take the limit n — > 4 
and extract the infinities in e R through the singularities of the T-function 
which appears in the series. The procedure leads to a renormalization mass 
counterterm of the type 
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As expected, the counterterm evaluated through the study of the lattice con- 
tinuum limit in n-dimensions contains a divergence in n for n — > 4 as well 
as a quadratic divergence for a — > 0. However, multiplying eq. (31) by a 2 , 
the second order pole in a disappears completely leaving the equation, now 
expressed in lattice units, finite (for n^4). 

The lattice derivation of the mass conterterm can be also performed directly in 
four-dimensions. Due to the important differences existing between the expan- 
sions in eq. (27) and eq. (28), a few changes take now place in the investigation, 



2 The explicit presence of the mass in the argument of the propagator coefficients 
has been included, in this case, simply to recall that it is m R rather than the bare 
m which now figures in the calculation. 
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though. The isolation of the infinities is no longer accomplished through the 
extraction of the poles in the T-functions, but resides ultimately in the obser- 
vation that the divergences in the lattice spacing mirror exactly the analogue 
divergent behaviours in e R . Operationally, the latter remark produces a four- 
dimensional lattice counterterm of the type 
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We notice that, in agreement with standard results [10], in this case both 
a quadratic and a logarithmic counterterm is needed at every order in the 
calculation. Nevertheless, in terms of lattice units, eq.(32) is, again, finite. 



4 Conclusions 



In the present work, we derived an asymptotic expansion for the modified 
Bessel function Y v {y 2 ($) as v — > oo. The expansion obtained was of vital im- 
portance to analytically evaluate the continuum expansion of both the lattice 
scalar propagator in a n-dimensional configuration space and its related gen- 
eralised quantity A L (x;n;g). The study of the small m(:r 2 ) ^-behaviour of 
A L (x; n; q) in the limit a — > was shown to involve only the standard series 
expansion of modified Bessel functions of either real (p rc ) or integer (p in ) or- 
der. The uniformity of the limits m — > and p rc — > p in was observed and the 
Liischer and Weisz expansion for the massless propagator in four-dimensions 
recovered. The result obtained in [4] for the mass-dependence was also repro- 
duced for integral values of p. Finally, as a perturbative application of the 
results obtained, the one-loop mass counterterm in A0 4 lattice theory was 
evaluated for both the n— and 4— dimensional case. 



12 



Acknowledgements 

Beatrice Paladini wishes to thank both the European Human and Capital 
Mobility Program and Hitachi Dublin Laboratory for their financial support 
towards the completion of this work. 

References 

[I] J. Collins, Renormalization (Cambridge University Press, 1984). 

[2] D. Z. Freedman, K. Johnson and J. I. Latorre, Nucl. Phys. B 371 (1992) 353-414. 

[3] M. Liischer and P. Weisz, Nucl. Phys. B 445 (1995) 429-450. 

[4] G. Burgio, S. Caracciolo, A. Pelissetto, Nucl. Phys. B 478 (1996) 687-722. 

[5] G. N. Watson, A treatise on the theory of Bessel functions, Second Edition 
(Cambridge University Press, 1962). 

[6] I. S. Gradshteyn and I. M. Ryzhik, Tables of integrals, series and products, 
Fourth Edition (Academic Press, 1965). 

[7] C. M. Bender and S. A. Orszag, Advanced mathematical methods for scientists 
and engineers (McGraw-Hill, 1978). 

[8] M. Creutz, Quarks, gluons and lattices (Cambridge University Press, 1983). 

[9] S. Wolfram, Mathematica (Addison- Wesley, 1988). 

[10] I. Montvay and G. Munster, Quantum fields on a lattice (Cambridge 
Monographs on Mathematical Physics, Cambridge University Press, 1994). 

[II] Z. X. Wang and D. R. Guo, Special functions (World Scientific, 1989). 

[12] A. F. Nikiforov and V. B. Uvarov, Special functions of mathematical physics 
(Birkhauser, 1988). 



13 



